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Particles held by springs in a linear shear flow exhibit oscillatory motion 
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The dynamics of small spheres, which are held by linear springs in a low Reynolds number shear 
flow at neighboring locations is investigated. The flow elongates the beads and the interplay of 
the shear gradient with the nonlinear behavior of the hydrodynamic interaction among the spheres 
causes in a large range of parameters a bifurcation to a surprising oscillatory bead motion. The 
parameter ranges, wherein this bifurcation is either super- or subcritical, are determined. 
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Introduction.- Studies about the motion of bacteria 
and flagella in a fluid, about the dynamics of blood 
cells and other small suspended objects such as poly- 
mers in simple flows are currently of central interest 
and count to one of the m ajor issues of microfluidics 
Qa!aiI0SIGl[!lE[l3 According to 
the short spatial scales involved in these cases the fluid 
motion surrounding the particles can be described in the 
small Reynolds numbers limit [j], 0, La- L| , by the linear 
Stokes equation [l|. 

The hydrodynamic interaction (HI) between neighbor- 
ing blood cells, swimming bacteria, different segments of 
a polymer or between neighboring polymers is of nonlin- 
ear nature even in the limit of low Reynolds numbers [1'1| . 
Moreover, this nonlinear behavior may cause dynamical 
effects, such as the periodic motion of small sedimenting 
spheres [l5lll6| | or the synchronization effects between ro- 
tating strings and between cilia [171 [18] , or it may cause 
a hydrodynamic coupling of particles in optical vortices 
19]. The HI may also amplify thermal fluctuations of 
polymers |2^ to mention only a few examples. 

Free single polymers show already a complex dynam- 
ical behavior in shear flow 0, an d the hydrodynamic 
interaction between many polymers leads to the so-called 
elastic turbulence [3, la] . Furthermore polymers fixed at 
one end in a plug flow are also a major issue [2(J, |2J| 
where one finds in this case significant hydrodynamic in- 
teraction effects both for the static as well as dynamic 
properties of the tethered polymers |2l], 122. l23j . Teth- 
ered polymers in shear flow have also been studied but 
so far only a single polymer fixed with one end at a wall 
was considered |2J, |25|, l'2oi |27| . Recently, investigations 
have been started in order to analyze the behavior of 
several flexible polymers fixed with their ends at the top 
of neighboring pillars |28| and exposed to a linear shear 
flow. So the interesting question arises quite naturally: 
what is the dynamics of neighboring tethered polymers 
in shear flow and which role plays the hydrodynamic in- 
teraction? 

We mimic a situation of interacting tethered polymers 
by spheres anchored by springs and neglect in a first ap- 



proach thermal fluctuations. To the best of our knowl- 
edge this is the first example where an oscillatory motion 
of bound particles in low Reynolds number flow has its 
origin in the hydrodynamic interaction. 
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FIG. 1: The three particle configuration in shear flow is 
shown. The corners of the triangle in the y — z-plane mark the 
minima of the harmonic potentials of the beads. The flow in- 
duced sphere displacements are indicated by the dashed lines. 
h measures the shift of the lower side of the triangle from the 
center of the shear flow, u(0) = 0. 

Model.- Three beads are fixed in a shear flow by linear 
springs with a spring constant k as shown in Fig.^ The 
locations of the minima of the corresponding harmonic 
potentials, Ri (i = 1,2,3), build, if not stated otherwise, 
an equilateral triangle of sidclcngth 6 and height H = 
b/\/2 with the upper corner at Ri = (0, h + H, 0) and 
the two lower corners at R2,3 = (0, /i, ±6/2). The bead- 
springs are elongated by a linear shear flow 



uo(y) = (72/,0,0) 



(1) 



and the actual bead positions r^ (i — 1,2,3) are deter- 
mined by the equations for the bead velocities 



uofo) - -fj 



J2(n,(v lJ )-n RP (r lJ )kv 3 )(2) 
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with fj = i\; — Rj. The first term describes the linear 
shear flow given by Eq. (|TJ and the second contribution is 
the ratio between the spring force and the Stokes friction 
C = Qnr]a that a single fixed particle of effective bead 
radius a experiences in a flow with shear viscosity rj. Each 




FIG. 2: The figure shows the stationary vertical displacements 
2/3 = 2/2 0I the lower beads (dashed line) and 2/1 of the upper 
bead (dash-dotted line) as well as the shift 23 = — z\ in z- 
direction (solid line) as a function of the vertical shift h/a 
and for the dimensionless shear rate T7 = 2.6. 



fixed particle causes a perturbation of the shear flow at 
the location of the other beads and vice versa. This so 
called hydrodynamic interaction (HI) is described for a 
Stokes flow by the Rotne Prager tensor [lj| 
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which describes together with the harmonic forces kij 
the fourth contribution in Eq. (J2J). ly = Sij is the unity 
matrix. The shear flow induces sphere rotations, which 
alter the flow field and therefore its action onto other 
spheres, as described by the third term in Eq. J2J 29] 
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where Ey = %{5i x 5j y + Si y Sj X ) if the particle can rotate 
freely and Ejj = "/5i y dj X if an external torque prevents 
the rotation. With the relaxation time r = C,Jk and the 
effective bead radius a one may rescale time t — + rt' , 
space r — > ar' and the shear rate 7' — > T7 and the results 
in this work are most conveniently presented in terms 
of these dimensionless units, as for instance the sphere 
displacement r d;i = r^ - R^ = (xd-i, Vd;i, Zd-,i) ■ 

Stationary displacement of the spheres.- The stationary 
solutions of the nonlinear equations @, i.e. r.j = 0, with 



the displacements r?.- := r' 



R' 



b d;i 
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-d;i. 



determined numerically by a Newton algorithm. r d 



are 
of a 

single bead increases according to the Stokes drag force 
F = Qitriau and the linear spring force linearly with the 
flow velocity Uo. However, by virtue of the nonlinear 
nature of the HI between the beads the elongation of the 
linear springs change nonlinearly as a function of the flow 
velocity, which itself varies linearly with the height h/a, 
as depicted in Fig. for rj = 2.6. 

For h = the flow velocity vanishes at the positions 
of bead 2 and 3 and with a finite shear gradient 7 only 
the upper bead 1 is displaced. The flow perturbation 



caused by bead 1 shifts beads 2 and 3 slightly downward 
and pushes both away in z— direction with z^. 3 — ~ z d-2- 
For finite h bead 2 and 3 are exposed to a finite veloc- 
ity Uo(r2 j 3) and excite also flow perturbations pointing 
both downward at bead 1. This involves j/°-i to become 
negative as well and since both perturbations act down- 
wards one has jy^ < y®. 2 at intermediate values of h, 
cf. Fig. 13 In z-direction both disturbances compensate 
each other, so that zj^ = is left unchanged. According 
to this stronger displacement y^. x at intermediate val- 
ues of h the relative distance between the upper and the 
two lower beads is reduced and therefore the flow pertur- 
bations caused by bead 1 are enhanced, and so are the 
values of |y^. 3 | and z^. 3 = —z®. 2 as a function of h. 

It is very surprising that all these displacements reach 
extrema, as shown in Fig.|2]and become smaller again for 
large values of h. An explanation of this behavior may be 
offered by inspecting the bead positions at large values of 
h/a. In this case the triangle built by the bead positions 
is again nearly parallel to the y — z plane and accordingly 
the flow perturbations and the related forces caused at 
the neighboring beads are nearly vanishing compared to 
their external force. In this limit, however, the height 
H of the triangle is smaller and the distance between the 
beads 2 and 3 is larger than for a vanishing fluid velocity. 
The latter behavior is a consequence of a complex bal- 
ance between the spring forces and the nonlinear forces 
due to the flow disturbances. Correspondingly there is 
hitherto no simple qualitative picture for both, the de- 
formed triangle built by the beads and the decreasing 
behavior of the elongation beyond their extrema. 
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FIG. 3: The largest real part a(h) is given in a) as a function 
of the vertical shift h/a for different shear rates 7 and part b) 
shows the corresponding imaginary parts u). 

Threshold of the Hopf bifurcation.- Slightly beyond the 
extrema in Fig. |3 the stationary bead displacements be- 
come unstable and one finds by numerically integrating 
Eqs. J2J using a standard method a bifurcation to oscil- 
latory bead motions. The threshold of this bifurcation 
may be determined by a linear stability analysis of the 
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stationary elongation r^ with respect to small pertur- 
bations Sri(t). Using the ansatz r^ = r" + 5ri(t), a lin- 
earization of Eqs. © leads to a set of 9 linear differential 
equations with constant coefficients 

Y = £(r°)Y with Y(t) = (6r u 6v 2 ,5r 3 ) , (5) 

governing the linear dynamics of the perturbations Sri (t) . 
Eq. J5J is solved by Y = exp(cr£' ± iwt')Yo, which trans- 
forms Eq. JSJ into an eigenvalue problem. The eigenvalue 
with the largest real part cr(h) has also a finite imagi- 
nary part lo and is positive within a finite range of h/a 
as shown in in Fig. [3] for three different values of the di- 
mensionless shear rate rj. The whole range of a positive 
cr(h) in the 7-7 — h/a plane is given by the shaded range 
in Fig. 0] The occurrence of oscillatory bead motion is 
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FIG. 4: a.) In the grey range the stationary bead elongation 
are unstable with respect to a supercritical Hopf bifurcation 
along the solid border line and to a subcritical one along the 
dashed line. Within the striped region the Hopf bifurcation 
is hysteretic. b.) The oscillation amplitude of bead 1 in x- 
direction is given as a function of h for 7-7 = 2.6. The dotted 
line marks the upper threshold of the Hopf bifurcation. 

rather robust with respect to changes of the anchor points 
of the linear springs which adhere the beads. We have 
tested this by changing the anchor point of bead I and 
2 in all three spatial directions. With such modifications 
the three anchor points build either no equilateral trian- 
gle or the equilateral triangle is inclined and not anymore 
perpendicular to the flow direction. The major trends in- 
clude the following ones. Bringing anchor points closer 
together enhances the hydrodynamic interaction which 
favors the Hopf bifurcation in a larger parameter range 
and it then takes also place at smaller shear rates and h. 
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FIG. 5: a.) The time dependence of the deviation r e; i from 
the center of mass for bead 1 (solid line) , 2 (dashed line) and 
3 (dash-dotted line) for 7-7 = 2.6 and h/a = 35.5. b.) The 
corresponding bead oscillations in the x — z plane. 



Nonlinear behavior of the bead oscillations.- A typi- 
cal example for the three dimensional oscillatory motion 
of the beads is given by a projection onto the x-axis 
in Fig. E^). Here the deviations r e .i = (x e - t i,y e - t i, z e -i) 
from the center of mass of the stationary solutions r cm = 
J2i=i 2 3 r ?/3 are displayed. Two characteristic features 
can be recognized. Bead 2 and 3 oscillate with a phase 
shift of 7r and bead I oscillates along the x-direction with 
twice the frequency of the other two beads. The double 
frequency of bead 1 is an effect of the projection onto the 
x-axis, as can be seen from the phase portrait in Fig.Et). 
Similar phase portraits can be obtained in the x — y and 
the y — z plane as well. Bead 1 performs a three dimen- 
sional motion and accordingly beads 2 and 3 are pushed 
away by a phase shift of tt as indicated in Fig.^Jr). 

The Hopf bifurcation is supercritical along the solid 
line bounding the grey range in Fig. 0Jl) . It is subcritical 
along the dashed one and the range of hysteresis is in- 
dicated by the striped range. The oscillation amplitude 
<5rf of bead 1 is shown in Fig. 0Jd) as a function of h/a 
at the shear rate T7 = 2.6. It indicates the supercritical 
behavior at the lower threshold and the hysteresis at the 
upper one. Close to the supercritical Hopf bifurcation 
the oscillations are harmonic. Further away from this 
threshold and in the parameter range with hysteresis in 
Fig-EJi) the periodic motion becomes rather anharmonic. 

Conclusions.- We found in this work a Hopf bifurcation 
of three bounded spheres in a low Reynolds number linear 
shear flow, which is induced by the interplay of the non- 
linear behavior of hydrodynamic interaction between the 
spheres and the shear gradient. To the best of our knowl- 
edge it is the first description of oscillations of bounded 



and hydrodynamically interacting particles in a Stokes 
flow. Most of the results are obtained for three beads 
anchored by linear springs at the corner of an equilateral 
triangle which is perpendicularly oriented with respect 
to the flow direction. The phenomenon is very robust 
against various variations of the anchor points. For two 
beads we did not find oscillations. 

Our results may also guide investigations on hydrody- 
namically interacting polymers fixed at small spheres and 
held by laser tweezers or anchored at boundaries in shear 
flow as well as for polymers that are fixed in shear flow 
close to boundaries at the top of pillars [28| . It is also an 
interesting question to be addressed, whether a rece ntly 
discussed cyclic motion for grafted polymers J2J, [2a, |22| 
is related to the Hopf bifurcation discussed here. 

We expect that several modifications of our model fa- 



vor oscillatory motion too. For instance nonlinear spring 
constants (which may mimic tethered polymers), differ- 
ent spring constants in different directions or an exposi- 
tion of the three beads to a Poiscuillc flow with its spa- 
tially dependent shear rate. The effects of these and other 
extensions are the subject of forthcoming work. 
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